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ABSTRACT 

Context. The detection of exolife is one of the goals of very ambitious future space missions or extremely large ground-based telescopes that 
aim to take direct images of Earth-like planets. While associations of simple molecules present in the planet's atmosphere (O2, O3, CO2 etc.) 
^ have been identified as possible global biomarkers, we analyse here the detectability of vegetation on a global scale on Earth's surface. 

, Aims. Considering its specific reflectance spectrum showing a sharp edge around 700 nm, vegetation can be considered as a potential global 
• • biomarker This work, based on observational data, aims to characterise and to quantify this signature in the disk-averaged Earth's spectrum. 
. 1^ Methods. Earthshine spectra have been used to test the detectability of the "Vegetation Red Edge" (VRE) in the Earth spectrum. We obtained 
reflectance spectra from near UV (320 nm) to near IR (1020 nm) for different Earth phases (continents or oceans seen from the Moon) with 
» I EMMI on the NTT at ESO/La Silla, Chile. We accurately correct the sky background and take into account the phase-dependent colour of the 
Moon. VRE measurements require a correction of the ozone Chappuis absorption band and Rayleigh plus aerosol scattering. 
Results. The near-UV spectrum shows a dark Earth below 350 nm due to the ozone absorption. The Vegetation Red Edge is observed when 
forests are present (4.0% for Africa and Europe), and is lower when clouds and oceans are mainly visible (1.3% for the Pacific Ocean). Errors 
are typically ±0.5, and ±1.5 in the worst case. We discuss the different sources of errors and bias and suggest possible improvements. 
Conclusions. We showed that measuring the VRE or an analog on an Earth-like planet remains very difficult (photometric relative accuracy of 
1% or better). It remains a small feature compared to atmospheric absorption lines. A direct monitoring from space of the global (disk-averaged) 
Earth's spectrum would provide the best VRE follow-up. 

Key words. Earth - Moon - Ultraviolet: solar system - Astrobiology - Techniques: spectroscopic - Vegetation red edge - Biomarker - 
Earthshine - Exolife 



1. Introduction 

Future space missions like Darwin/TPF will provide the first 
images and low-resolution spectra of Earth-like extrasolar 
planets, with the final goal to extract a "biosignature" from the 
spectra. The TPF-C NASA mission will observe exo-terrestrial 
planets around the nearby stars with a large off'-axis telescope 
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equipped with a visible coronagraphic device. To prove the 
interest of the visible part of the spectrum in spite of the 
huge contrast (10'°) compared to thermal infrared (10^), we 
study the presence of biomarker features. The spectrum of 
the light reflected by a planet, when normalized to the parent 
star spectrum, gives the planet reflectance spectrum, reveal- 
ing planet atmospheric and ground colour, when the latter is 
visible through a partially transparent atmosphere. The planet 
will remain unresolved and its spectrum will be spatially in- 
tegrated (i.e. disk-averaged) at the observed orbital phase of 
the planet. In this work, we estimate how the Earth's spectrum 
would look like from ~ 10 pc. This can be done - in prin- 
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ciple at least - by integrating spatially-resolved spectra from 
low-orbit satellites. It can also be done from a space probe 
traveling into the Solar System and looking back the Earth, as 
Mars Express did in 2003' . An alternative method to obtain the 
Earth's disk-averaged spectrum consists of taking a spectrum of 
the Moon Earthshine, i.e. Earth light backscattered by the non- 
sunlit Moon. A spectrum of the Moon Earthshine directly gives 
a disk-averaged spectrum of the Earth at a given phase as seen 
from the Moon (since the Moon surface roughness "washes 
out" any spatial information on the Earth's colour). 

Biosignatures are of two types: out-of-equilibrium 
molecules in the planet atmosphere (like oxygen and 
ozone) or ground colours characteristic of biological com- 
plex molecules (like pigments in vegetation). The visible 
and near infra-red Earthshine spectra published to date 
(Arnold et al. 2002 Woolf et al. 20Q2| [Seager et al. 2005] 
[Mon tanes-Rodriguez et al. 2005 1 clearly show the atmo- 
spheric signatures and, at least, tentative signs of ground 
vegetation which thus appears as an interesting potential 
global biomarker Vegetation indeed has a high reflectivity 
in the near-IR, higher than in the visible by a factor of ^ 5 
(Clark 1999 1. This produces a sharp edge around x 700nm, the 
so-called Vegetation Red Edge (VRE). An Earth disk-averaged 
reflectance spectrum can thus rise by a significant fraction 
around this wavelength if vegetation is in view from the Moon 
when the Earthshine is observed. Arcichovsky suggested in 
1912 to look for chlorophyll absorption in the Earthshine 
spectrum, to calibrate chlorophyll in the spectrum of other 
planets (Arcichovsky 1912 1. This approach re-emerged again 
only in 1999 with preliminaryt tests at ESO by several among 
us (RE, J.S. and D.B.) and at OHP (Arnold et al. 2002i. The 
position of the VRE may not lie at 725 nm on an exo-Earth. Its 
possible confusion with mineral spectral features is discussed 
by Schneider (2004) and Seager et al. (2005). 

The red side [600:1000 nm] of the Earth reflectance spec- 
trum shows the presence of O2 and H2O absorption bands and 
of the VRE, while the blue side [320:600 nm] clearly shows 
the Huggins and Chappuis ozone (O3) absorption bands. The 
higher reflectance in the blue shows that our planet is blue 
due to Rayleigh scattering in the atmosphere, as detected by 
Tikhoff' (1914) and Very (1915), and confirmed later with ac- 
curate Earthshine photometry by Danjon (1936). 

The present paper presents new observations and em- 
phasizes the difficulties of data reduction and describes 
how it has been improved with respect to our first work 
JArnold et al. 200 2 1. In Section ISHI we explain how we adapt 
Qiu's method ( |Qiu et al. 2003t to carefully subtract the sky 
background from Earthshine spectra. Section l331 describes how 
we correct the phase-dependent colour eff'ect of the Moon, fol- 
lowing a new processing of Lane & Irvine (1973) data. This 
correction was not done in Arnold et al. (2002). Section 0] 
presents the first - as fas as we know - near-UV integrated 
Earth spectra obtained from Earthshine observations that reveal 
a dark Earth below 350 nm due to the strong O3 absorption 
(ozone Huggins absorption bands). While this work presents 
smaller VRE than Arnold et al. (2002) and Woolf et al. (2002), 

' http://www.esa.int/esaCP/Pr. 44_ 2003_ p_ EN.html 



it nevertheless shows different values depending whether ocean 
or land are seen from the Moon at the epoch of observations. 

2. Observations 

The observations were taken at the NTT/La Silla telescope 
(3.5m) on July 24th and September 18th 2004 for the descend- 
ing Moon and May 3 1 st and June 2nd 2005 for the ascending 
Moon (TableQ. The Earthshine spectrum for 05-31-2005 was 
taken just before the last quarter of the Moon and has a very low 
contrast compared to the sky background. The spectra were ob- 
tained with the EMMI spectrograph in the medium dispersion 
mode in the blue (BLMD mode) and in low dispersion in the 
red (RJLD mode), enabling us to record spectra from 320 to 
1020 nm, with a 20 nm gap around 520 nm. The spectral res- 
olution is R=450 in the blue and R=250 in the red. To record 
both Earthshine (hereafter ES) and sky background spectra si- 
multaneously, EMMI's long slit was oriented East- West on the 
lunar limb ( Arnold et al. 2002 1. ES exposures are bracketed by 
at least two exposures of the bright Moonshine (hereafter MS) 
with the slit oriented North-South. The length of the slit (6- 
arcmin in blue and 8-arcmin in red modes) allows us to sample 
the Moon spectrum over a large lunar region, giving a correct 
mean of the Moon spectrum. MS spectra are recorded through 
a neutral density filter in the blue arm and with a diaphragm in 
the red arm (unfortunately the diaphragm could not be placed 
exactly in a pupil plane, resulting in strong vignetting - no neu- 
tral density filter was available in the red). 

3. Spectra processing method 

Data reduction is done with dedicated IDL^'^routines. All im- 
ages are processed for cosmic rays, bias, dark, flat corrections 
(dome-flat) and distortion. Some of the MS blue spectra have 
been obtained through a neutral density filter (ND=0.3) and 
these files have been corrected from its not-perfectly-flat spec- 
tral profile. The Earthshine spectrum needs to be corrected 
for the bright sky background (comparable to the Earthshine 
flux), the Earth atmosphere transmittance and phase-dependent 
colour effects of the Moon reflectance. 

3.1. Sky background subtraction 

A fraction of the sunlit Moon flux is scattered by the atmo- 
sphere and produces a bright background around the Moon, 
decreasing with the distance to the sunlit crescent. Both 
sky and Earthshine must be corrected from this "pollution". 
We consider that this noise decreases linearly with distance 
(Qiu et al. 2003). The sky background spectrum recorded near 
the ES spectrum can thus be linearly extrapolated to find its 
value over the Moon surface. We first make a linear extrap- 
olation at each wavelength to construct what we call a syn- 
thetic sky background to be subtracted from the Earthshine 
(Fig. [0. The figure is in agreement with the linear assump- 
tion mentioned above. If we plot the slopes obtained at each 
wavelength for the linear extrapolation, we notice strong varia- 
tions between adjacent wavelengths. These variations are pro- 
portional to the flux at each wavelength, but are also due to 
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Date 


Hour 


Exposure Time 


Exposure Time 


MS phase angle 




(UT) 


(blue) 


(red) 




05/24/2004 


23hl6(blue) 


2xl80sec. 




-117° 


09/18/2004 


Oh 15 (red) 




2xl00sec. 


-137° 


05/31/2005 


9h30(red) and lOhOl(blue) 


2x250sec. 


2xl20sec. 


101° 


06/02/2005 


9h25(red) and 9h50(blue) 


2x250sec. 


2xl20sec. 


126° 



22 
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Fig. 1. Flux at a given wavelength along the slit placed per- 
pendicular to the edge of the Moon to record both Earthshine 
and sky background simultaneously. The linear fit based on 
the sky background is extrapolated through the Moon surface 
to estimate the synthetic sky background to be removed from 
the Earthshine data. The transition zone between ES and back- 
ground (pixel 200 to 280) is not used in the data processing. 




800 
ambda in nm 



Fig. 2. The slopes obtained at each wavelength from the sky 
background (Fig. [0 are normalized to the mean flux of the 
recorded sky background (dashed line) and smoothed (solid 
line). 



3.2. The Earth reflectance ER{A) 

The Earth reflectance ER is given by the ratio of the 
Earthshine spectrum divided by the sunlit Moon spectrum 
dArnold et al. 2002t 



MS(A) 



(1) 



Note that the ES spectrum contains the signature of roughly 
three paths through the Earth atmosphere: Two paths results 
from direct Sun light that is afterward scattered toward the 
Moon and a thkd one from the light coming back from the 
ES. MS results from only one path. Therefore the division in 
Eq[nexactly corrects one path through the atmosphere, which 
is identical for both ES and MS spectra. Thus ER contains the 
absorption signature of roughly two paths through the atmo- 
sphere (see lArnold et al. 2002l for details). 

Eq^assumes that ES(/l) and MS(/t) are recorded simulta- 
neously, i.e. at the same air-mass. In practice, each ES exposure 
is bracketed by two MS exposures. To estimate the MS spec- 
trum at the epoch of the ES exposure, the MS is obtained from 
the average of the two MS spectra, weighted proportionally to 
the time elapsed before and after the ES exposure. For the blue 
spectra of the MS, we also need to correct for the effect of the 
neutral density filter that does not have the same absorption 
at each wavelength. The r{A) is a geometrical chromatic factor 
that takes into account the geometrical positions of the Sun, the 
Earth and the Moon. It is described in the next section. 

A photometrically calibrated Earth reflectance spectrum 
can be obtained in principle by following the procedure used 
for broad-band albedo measurements ( |Qiu et al. 20 03). This re- 
quires us to measure spectra or brightness only in calibrated 
areas to properly take into account the Moon phase function. 
Our data do not allow us to calibrate all fluxes, so our Earth 
reflectance spectra are not absolutely calibrated. 



the lower S/N in the absorption lines. To reduce the noise on 
the strongest telluric lines, the slopes of sky background are 
normalized to the mean flux at each wavelength and smoothed 
(boxcar average of 4 to 20 pixels width depending on the data. 
Fig. |2j- The final synthetic sky background spectrum is ob- 
tained by multiplying the previous smoothed spectrum by the 
sky background mean flux at each wavelength. Once the back- 
ground is corrected, each image is binned into an improved S/N 
ratio spectrum. 



3.3. The Moon phase colour effect 

Lane & Irvine (1973) described the chromatic dependence of 
the integrated Moonshine versus the lunar phase. Their narrow 
band photometry showed a ^30% excess of red light at 1050 
nm with respect to the flux at 350 nm at phases 120°. Lane & 
Irvine plotted MS fluxes versus phase angles for different wave- 
lengths. They identified a linear domain for low phase angles 
(Iphasel < 40°) and applied a cubic fit for |phase| > 40°. Their 
cubic fit is strongly constrained by a very few data points at 
high phase angles (around 120°), which probably biases their 
fit. They also observed a difference of brightness between the 
waxing and waning Moon observed at exactly opposite phase 
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60 80 
phase angle 

Fig. 3. Lane & Irvine data at 1063.5 nm: stars and crosses cor- 
respond to data from 1964 and 1965, respectively. The solid 
line is the original Lane & Irvine linear (<40° ) and cubic (>40° 
) fit based on the 1965 set of data only. The dash line is the sec- 
ond order fit based on all data and used in this work to correct 
the phase-dependent colour of the Moon. 



600 700 
arrbda In nm 



Fig. 4. Data extrapolated at a Moon phase of -137°. The crosses 
are determined with a linear-cubic fit and the star with our sec- 
ond order fit. The data are normalised at 1020 nm for easier 
comparison. A second order fit is used to define the colour cor- 
rection (plain line for the star and dotted line for the crosses). 



angles. This is due to the different proportion of visible craters 
and Mares ( |Rougier 1933| IRussell l"916 ) with respect to the 
Moon phase. But Lane & Irvine report that the difference they 
measured seems too large. They processed only their second set 
of data, with the majority of measurements at positive phase an- 
gles. This slightly biased their results toward positive phase an- 
gles, but their data were nevertheless consistent with Rougier's 
data (1933), the most accurate at that time. 

Applying the Lane & Irvine correction to our data leads to 
a positive slope of the ER spectrum in the [400:520 nm] do- 
main. This is inconsistent with the dominant Rayleigh scatter- 
ing which gives a negative slope in the blue for a cloud-free 
Earth, or, at worst, a zero slope for a white Earth fully cov- 
ered by clouds. Paillet & Selsis (2005, private communication) 
confirmed the permanent negative slope in the blue with simu- 
lated Earth reflectance spectra obtained with a radiative transfer 
code. Therefore, the Lane & Irvine original fit leads to a clear 
over-correction of the colour phase eff'ect. 

We therefore reconsidered all the Lane & Irvine data and fit 
all points with a second-order curve over the full range of phase 
angles, without considering two distinct regimes (i.e. linear + 
cubic). An example of the fit is shown in Fig.|3] We processed 
each wavelength and extrapolated the correction at the phase 
angles of our observations. Another second-order interpolation 
through the nine obtained values (the nine Lane & Irvine pho- 
tometric filters) gives the r{A) factor needed to correct the re- 
flectance spectra given by Eq.[n(Fig.|4}. 

We have to observe two regions of the Moon, the ES and 
the MS. While MS is observed at high phase angles (i.e. a nar- 
row crescent), ES is observed at 0° phase angle because the 
Moon source of light and the observer are both on Earth. This 
difference in phase angles induces a colour difference between 
ES and MS that must be taken into account in the data reduc- 
tion. This is done with the r{A) factor (Eq.[l) which allows one 
to calculate the MS spectrum at 0° phase angle, i.e. the same as 
forES. 

We also compared each individual ES spectrum line along 
the slit going from the edge to the center of the Moon and de- 
tected a smooth chromatic variation along the slit, resulting in 



a redder Moon at the center than at the edge of the lunar disk. 
This residual chromatic variation is due to i) spectral changes 
due to the soil nature (crater versus Mare), and ii) a chromatic 
dependence of the spectra along the spectrograph slit, due to 
the Bidirectional Reflectance Distribution Function (BRDF) of 
the soils sampled along the slit. Most of the phase-dependent 
colour effect is canceled by the Lane & Irvine correction, but it 
provides only a mean correction based on observations of the 
integrated illuminated Moon at a given phase, while here we 
speak about differences from point to point on the lunar surface. 
In our spectra, the maximum increase in the red between 320 
and 1020 nm is of the order of 3%. The measurement of this 
residual colour effect however remains too uncertain. It seems 
smooth and minor so we have neglected it in our data reduction, 
although we estimate an induced bias of +0. 15 on the VRE. 

3.4. The VRE measurement method 

To detect the vegetation signature, i.e. a signature of the 
ground, in principle we need to extract the Earth surface re- 
flectance from the ER, which contains the signatures of the 
Earth atmosphere. To remove the atmospheric signatures, ER 
can be divided by the atmosphere transmittance AT(/l) that can 
be obtained from the ratio of two spectra of a calibration star (or 
MS) taken at two different air-masses (Arnold et al. 2002jl. It 
can also be built from spectral databases and adjusted to fit the 
observed spectrum IWoolf et al. 2002t . In this work, we only 
correct the ER spectrum for the Oj, absorption because, among 
the main atmospheric absorption lines (O2, and H2O), only 
the broad Chappuis band extends up to w 700 nm and thus can 
potentially induce a bias in the red part of the spectrum where 
the VRE is measured (Fig.|5}. 

The ER(/l) spectra are fitted with laboratory measured 
spectra ( |Voigt et al. 200 1^ convolved to reproduce EMMI's 
spectral resolution. Removing of the Chappuis band re- 
quires great care: we observed that the band can indeed 
easily be slightly over-corrected, leading to an apparently 
smooth Rayleigh scattering but an underestimated VRE. In the 
[520:670 nm] domain, the ER spectrum is dominated by the 
absorption of Oj,, Rayleigh and aerosol scattering. To remove 
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Fig. 5. Spectra of the main atmospheric absorbing components 
(shifted vertically for clarity). Vertical lines indicate the two 
spectral bands used to calculate the VRE, [667:682 nm] and 
[745:752 nm], and show that only O3 intersects the domain 
where the VRE is measured. 
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Fig. 6. ER(/l) spectrum (dotted line) and corrected for the O3 
Chappuis absorption band (solid line). The smooth solid line is 
a fitted Rayleigh-aerosol scattering model. Vertical lines define 
the two spectral bands used to calculate the VRE. 



the O3 Chappuis band, we fit a Rayleigh-aerosol function based 
on continuum parts of the spectrum for different concentration 
of O3. We then fit a second-order polynomial on the root mean 
squares of the difference between the fit and the spectrum for 
each Oj, concentration to find the best fit, i.e. the minimum root 
mean square and its corresponding O3 concentration. The hor- 
izontal RMS of this second fit gives an estimate of the ozone 
column. We define the Rayleigh-aerosol function as the sum of 
a Rayleigh scattering function a + bA^* plus an aerosol scat- 
tering contribution c/l"'^. To reduce the number of parameters, 
we consider that the aerosol and the Rayleigh contribution are 
equal at A,, = 700 nm JLena 1996b . so the scattering law is 

scat(^) = a-.b(l + -^). (2) 

Once the Chappuis band is corrected, the spectrum is again fit- 
ted with a Rayleigh-aerosol function in continuum domains and 
normalised to that function to obtain a spectrum where the veg- 
etation signature can be quantified, although the spectrum is 
not corrected for O2 and H2O. The VRE is measured on this 
corrected spectrum (Fig.|6} following Eq|3 

VRE = ^ (3) 

where and r/ are the mean reflectance in the [667:682 nm] 
and [745:752 nm] spectral bands. 

4. Results and discussion 

4.1. Qualitative analysis 

Figures |H1 and |51 show the obtained ER spectra. The spectra 
show different features: 

- The blue part (Fig|S|l displays a decrease above 350 nm due 
to the Rayleigh scattering. This indicates a partially clear 
atmosphere and shows that the integrated Earth could be a 
'pale blue dot' or a nearly white dot (2nd June 2005 data), 
given the amount of cloud. 

- At wavelengths below 350 nm (Fig|5), reflectance becomes 
much lower indicating that Earth is clearly a much fainter 
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Fig. 7. Resulting ER(/l) corrected for O3 Chappuis absorption, 
Rayleigh and aerosol scattering, thus ready for the VRE mea- 
surement. Vertical lines define the two spectral bands used to 
calculate the VRE. The plots have been shifted vertically for 
clarity. 

object in the near UV. This is due to the ozone strong ab- 
sorption in the UV (Huggins band) giving a high contrast 
between the visible and UV Earth. The three UV spectra 
are slightly different due to the different signal to noise ra- 
tio of the data. It could be possible to measure the integrated 
ozone column from the Huggins band ( Arnold et a l. 2002> . 
but we could not investigate this in more details due to the 
low signal to noise ratio is this spectral domain, cause by 
the low signal and residual solar spectral features that are 
not perfectly cancelled out in the division of Eq. ^ The 
wider but weaker band of O3 absorption (Chappuis band) 
peaks at about 600 nm. 
- Oxygen bands at 630, 690, and mainly at 760 nm 
are clearly visible. The oxygen dimer (02)2 produces a 
weak absorption band visible between 570 and 580 nm 
frinetti et al. 20061 1. but is superimposed by a weak ozone 
features from the Chappuis absorption band. The depth of 
the 760nm band is routinely used in the Earth remote sens- 
ing to measure the altitude of clouds by satellites imaging 
Earth's surface: the deeper the absorption, the lower the al- 
titude of the cloud. Thus in principle the depth of this band 
in integrated spectra gives an average of the clouds altitude. 
On an extra-solar planet, the band's relative depth variation 
with time would indicate a variation of cloud mean altitude. 
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Table 2. Vegetation Red Edge. The Ot, concentration is the 
amount that provides the best (least-square) correction of the 
Chappuis band. The raw VRE is compute on the reflectance 
spectra without any atmospheric correction. 



Date 2004/09/18 2005/05/31 2005/06/02 



VRE[%] 1.3+0.5 3.4+0.5 4.0±0.5 

O3 [DU]±3o- 1100±600 110±2 740±70 
rawVRE[%] 1.4 2.0 5.8 

In our data, the relative depth of O2 shows a larger variation 
than expected, for a not yet understood reason. 

- Water vapour bands are also clearly visible features, espe- 
cially in the redder part of the spectrum. 

- Vegetation produces a small increase in the continuum 
spectrum, barely visible in the 730 to 750 nm range com- 
pared to the bluer continuum, and requiring a more quanti- 
tative analysis described in the following section. 

Since O2 and Ot, are product and by-product of photosyn- 
thetic life on Earth, the observation of these two molecules in 
the spectrum of an extra-solar planet could be a potential indi- 
cator of life, in combination with other parameters (planet in 
the habitable zone, i.e. with liquid water). 

4.2. Quantitative analysis 

During the first observing run, both blue and red spectra were 
not recorded during the same night, but the aim was to record at 
least a red spectrum of the Earth reflectance showing the veg- 
etation (09/18/2004). The second run allowed us to record all 
spectra from 320 to 1020 nm in about 40 min. We measured a 
VRE between 1.3 and 4.0% (Table |2 Figure with an accu- 
racy of about +0.5 (absolute error on VRE). 

The many steps in the data reduction process required to 
obtain ER and VRE induce the uncertainties evaluated as fol- 
lows: 

- Small variations in the width of the continuum used for the 
determination of the VRE give variations of about +0. 1 (i.e. 
VRE = 4.0 + 0.1%). This comes from the noise in the spec- 
trum over the relatively narrow domains of continuum. 

- The fit of the Rayleigh-aerosol function is based on narrow 
continuum domains between 520 and 680 nm. Variations 
in the limits of these domains result in a small +0.05 on 
the VRE, but in the noisy spectra of 05-31-2005, the varia- 
tions can be up to +1.5. This shows the care needed in the 
acquisition and processing of the data. 

- Table|2]gives the amount of O3 giving the best (least square) 
Rayleigh-aerosol fit of the data. The mean terrestrial O3 
column density is known to be about 300 Dobson Units 
(DU)^. Considering that the reflectance ER results from 
two paths through the atmosphere and is integrated on the 

^ A Dobson unit (DU) represents the thickness of a layer of an 
atmospheric gas, here ozone, if it could be compressed to standard 
ground temperature and pressure. One Dobson Unit is 0.01 mm and 
the total ozone column is typically 300 DU. 



Solar Irradiance 




320 325 330 335 340 345 350 

ambda in nn^ 

Fig. 9. ER(/l) near-UV spectra (three lowest curves) showing 
the strong absorption of O3 ozone (Huggins band). The theo- 
retical spectrum for ozone is from Voigt et al. (2001). The up- 
per spectrum is the solar irradiance JArvesen et al. 1969t . The 
noise in the three near-UV ER spectra is due to the low signal 
in the UV and to imperfect correction of the solar irradiance 
by the division of Eq. [J The spectra are shifted vertically for 
clarity (-1-O, H-0.3, -hO.6, -i-l.l, +1.6). 

visible illuminated Earth (see Sect l3.2> . we expect to mea- 
sure an equivalent column of the order of 1200 DU. But 
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Fig. 8. Earth reflectance spectra. The blue spectra have been scaled to agree (least mean square) with the Rayleigh-aerosol fit of 
the red spectra. On the right, the Earth seen from the Moon at the epoch of the observation (cloudless model). 
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we find values from 1 10 to 1 100 DU (Tabled), suggesting 
that the Ot, fit based only on the narrow continuum parts of 
the red side of the ER spectrum is probably not sufficient. 
A better fit would have been obtained with a full spectrum 
from 400 to 800 nm to benefit from the full Chappuis band 
for the fit, but unfortunately it was not possible to overlap 
both blue and red spectra. The estimation of the Oi concen- 
tration results from the minimum root mean square between 
the fit and the data. The uncertainty in the ozone correction 
gives a variation of ±0.2 on the VRE. 

- The determination of the synthetic sky background is done 
by smoothing the variations of the background slope per 
ADU. The width of the smoothing box is determined by 
the S/N ratio in the ES spectrum. Choosing a different but 
still reasonable smoothing box produces variations of +0.1 
on the VRE, and if we do not smooth the values, the VRE 
varies by around +0.5. 

- The phase-dependant colour of the Moon is estimated from 
a relatively reduced and sparse number of data points from 
Lane <& Irvine (1973), but as far as we know, these are 
the only available at phases above 90°. Comparing differ- 
ent approaches for the extrapolation of the Lane & Irvine 
data to larger phase angles, the VRE remains unchanged 
within ±0.1. But if we do not take into account the phase- 
colour correction, all VRE values decrease by about 1.5, 
pointing out the relevance of the colour correction. The 
residual colour bias of 3% (end of Sect l3.3> induces a 
positive bias in the VRE of a; 0.15. To minimize the 
colour bias, better observations should always observe the 
same calibrated spots in Earthshine and sunlit Moon to ob- 
tain the exact phase-dependent colour and photometry of 
these areas, as Qiu did for accurate albedo measurements 
<Qiuetal. im^ . 

Simple models (Arnold et al. 20021 'Arnold et al. 2003 1 
predicted a VRE near 8-10% when Africa and Europe light the 
Moon. Our results, from 3 to 4%, are lower, but it is necessary 
to have more spectra to validate our treatment method or com- 
pare with simulations including the real cloud coverage. The 
result of 1.3% of VRE for the 09-18-2004, when the Pacific 
faced the Moon, suggests that there is probably a small posi- 
tive bias in our VRE estimation method. We also note a small 
positive slope in the red part of the red spectrum, between 700 
and 1000 nm, suggesting the signature of deserts (Sahara). If 
we correct the spectrum for this increase, all VRE estimations 
decrease, giving a VRE near 3% for the Earth showing Africa 
and Europe. This low VRE may be consistent with the satel- 
lite images that show large cloud coverage over rain forests. 
However, VRE over the Pacific Ocean versus Africa remains 
different thus allowing the detection of vegetation on Earth. 

As described in Section 13.31 the best way to correct for 
the phase-dependent colour effect of the Moon is to have cali- 
brated photometric data. When we record the Earthshine, we 
multiplied the reflectance spectra of the Earth by the Moon 
reflectance and the atmospheric transmittance. The Earthshine 
is always observed at 0° phase. To correct for the Moon re- 
flectance, we thus need a record of the same pattern at Full 
Moon, corrected from the atmosphere transmittance. To cor- 



rect for the atmosphere transmission during the observations, 
we need the spectrum of the Moonshine with a calibrated re- 
flectance for that particular phase. 

5. Conclusion 

We have presented spectra of the Earth reflectance from near- 
UV to near-IR (320 to 1020 nm) for four different nights. 
They show significant variations in Rayleigh scattering de- 
pending on the cloud cover (the Earth 'blue dot' can be al- 
most white). One of the spectra was taken with mainly light 
from the Pacific Ocean while the others included parts of or 
the whole Africa and Europe. The Vegetation Red Edge is ob- 
served when land with forests are present, with values between 
3 and 4% (± « 0.5%), but remains very low otherwise (^ 1%). 
A spectro-photometric accuracy better than 1 % is therefore re- 
quired for future ground based observations (with Extremely 
Large Telescopes) or space missions aiming to detect an analog 
to the terrestrial VRE. These instruments will all be equipped 
with an optimized coronagraph allowing starlight rejection bet- 
ter than 10^ to be able to provide the image of the planet itself. 

A survey over one year or more with monthly observa- 
tions would allow one to follow the VRE seasonal variations 
and undoubtedly improve our knowledge of the behavior of 
this biomarker. In addition, it would be desirable to have a di- 
rect measurement and monitoring from space of the global (i.e. 
disk-averaged) spectrum of the Earth. One solution would be 
to observe Moon Earthshine: from space, the Earth reflectance 
would still be obtained with EqQ but would avoid the source 
of noise from one path through the atmosphere for both ES 
and MS light. Other solutions would be to integrate spatially- 
resolved spectra from low-orbit satellites, or ideally, to have a 
satellite far enough from the Earth to see it entirely and able to 
easily obtain disk-averaged spectra of the Earth radiance. The 
Earth reflectance would then simply be obtained after a divi- 
sion by the solar radiance, avoiding all sources of noise and 
bias due to the Moon. 

Our work showed that measuring the VRE remains very 
difficult, and, although measurable, the VRE is a small feature 
when compared to O2, O3 and water vapor absorption lines. 
But O2 present in the Earth atmosphere is 100% of photosyn- 
thetic origin and its deep absorption line at 760 nm is, for the 
Earth, a signature of photosynthesis. This suggests that if O2 
is observed at 760 nm on an exo-Earth, it is relevant to look 
for a VRE-like (probably small) feature in the spectrum, keep- 
ing in mind that this signature can be significantly different 
(spectral shape and wavelength) on an exo-Earth than on our 
planet. Clearly a sharp feature may be much more detectable 
than a smoother signature. An exhaustive analysis of possible 
artifacts remains necessary to consolidate the work done on 
this subject CSchneider 200 4 Seager et al. 2005^ . Phylogenetic 
analysis of plants and research on alternative photosynthesis 
chemistries are also of prime interest to better constrain this 
plausible biomarker. 
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